Hysteresis and intermittency in a nano-bridge based suspended DC-SQUID 
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We study voltage response of nano-bridge based DC-SQUID fabricated on a Si3 N4 membrane. 
Such a configuration may help in reducing 1/f noise, which originates from substrate fluctuating 
defects. We find that the poor thermal coupling between the DC-SQUID and the substrate leads 
to strong hysteretic response of the SQUID, even though it is biased by an alternating current. In 
addition, when the DC-SQUID is biased near a threshold of spontaneous oscillations, the measured 
voltage has an intermittent pattern, which depends on the applied magnetic flux through the SQUID. 

PACS numbers: 



Superconducting quantum interference devices 
(SQUIDs) are key components in many applications 
One of the major research goals with SQUIDs nowadays 
is to reduce noises in the Josephson junctions (JJs) 
composing the SQUIDs. This issue was extensively 
studied in the past decades, where it was shown that one 
of the main sources of noise is the couplingof the JJs 
to the substrate, which induces 1/f noises [1J. Several 
studies have tried to solve this problem by suspending 
the JJs and thus decouple them from the substrate. 
So far these studies have only shown moderate success 
Another emerging technology, which attracts an 
increasing interest, is SQUIDs based on nano-bridge 
weak-links. Such SQUIDs can be made extremely small 
and have the potential to outperform conventional 
SQUIDs in terms of noise properties @, 0] • 

Our original research goal was to study noise properties 
of a nano-bridge based DC-SQUIDs suspended on a Sili- 
con Nitride (S13 N4) membrane. However, we discovered 
that non-equilibrium thermal processes, which emerge 
due to the poor coupling of the JJs to the substrate, 
play a dominant role in SQUIDs dynamics. When the 
SQUID is excited by alternating current strong hysteretic 
response of the DC-SQUID is observed, even though such 
excitation should, in principle, eliminate hysteretic be- 
havior. In addition, when the SQUID is biased near a 
threshold of spontaneous oscillations, the measured volt- 
age shows an intermittent pattern, which depends on the 
magnetic flux through the SQUID. To account for these 
results we extend the theoretical model in Ref. @ to 
include low frequency variations in the SQUID tempera- 
ture. Numerical simulations show qualitative agreement 
with the experimental data. 

A simplified circuit layout of a typical device is illus- 
trated in Fig. DJb). We build our devices on a high 
resistivity Silicon wafer, covered by a 100 nm-thick layer 
of Si3 N4. A small section of the wafer is etched from the 
back to produce a suspended Si3 N4 membrane, having 
size of 100 /1m 2 . The DC-SQUID is made of Niobium, 
having layer thickness of 80 nm and lateral dimensions of 



300K (b) 4.2K 

I fSilicon Wafer 1 



Lock-in -VW- 
Amplifier 
IMC) -VW* — 



lkQ 



Oscilloscope 
1MQAWK 



-Hlpf 



I 



si 



.A 



IjjLPFi ■ j_ 




* Current affiliation: Department of Physics, University of Basel. 



FIG. 1: (Color online) (a) Measurement setup, (b) A simpli- 
fied circuit layout for a DC-SQUID, (c) Electron micrograph 
of a DC-SQUID partially fabricated on top a Si3 N4 mem- 
brane, (c) Electron micrograph image of a nano-bridge. 



110 x 7/im 2 . The self-inductance of the SQUID is calcu- 
lated using FastHenry program Q to be L = lOOpH. The 
DC-SQUID is composed of two nano-bridge JJs (NBJJs), 
one NBJJ in each of its two arms. The dimensions of the 
bridges are 100 x 115 nm 2 , and their combined critical 
current is 1.2 mA. Most of the SQUID structure, includ- 
ing one of the NBJJs is fabricated on the Si3 N 4 mem- 
brane, as shown by the electron micrographs in subfig- 
uresQJc) andQJd). An on-chip stripline passes nearby the 
DC- SQUID, and is used to generate magnetic flux in the 
DC-SQUID loop. A DC bias line, which includes on-chip 
low-pass filters (LPFs), is connected to the DC-SQUID 
and is used for voltage measurements of the SQUID. Note 
that the DC-SQUID is embedded in a superconducting 
stripline resonator (not illustrated in Fig. [ljb)), but the 
resonator was not used in the experiments, and thus had 
a negligible influence on the experiments presented in this 



2 



E 60 (a) 

§20- co/27i=10kHz 



0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 
I [mA] 



< 
E 



0.8 1 





(0) 








.itic 






th 






•dec 






th 











10 



10 10 
co 12k [Hz] 



0.98 
-° 0.96 
jS 0.94 

0.92 





"(c) 


.inc 






th 






.dec 






th 











1 ffl , c M 



10 



FIG. 2: Hysteretic current-voltage measurements. Panel (a) 
shows DC-SQUID voltage traces measured for increasing and 
decreasing current excitation sweeps. Panels (b) and (c) show 
measurement and simulation results, respectively, of increas- 
ing and decreasing threshold currents as a function of the 
excitation frequency. 



paper. Further design considerations and fabrication de- 
tails can be found elsewhere Gil [111 . 



Our experiments are carried out using the setup de- 
picted in Fig. [TJa). We use a lock- in amplifier, which 
applies alternating current through the SQUID. The ex- 
citation frequencies can be up to 100 kHz. We measure 
the voltage across the DC-SQUID using the lock-in am- 
plifier and record its time domain dynamics using an os- 
cilloscope. In experiments showing intermittency we also 
apply DC magnetic flux through the SQUID. All mea- 
surements are carried out while the device is fully im- 
mersed in liquid Helium. 

Figure Uta) shows an example of the DC-SQUID hys- 
teretic response to an excitation current having alter- 
nating frequency of w x /27r = 10 kHz. This measure- 
ment shows the measured voltage across the SQUID 
when the current amplitude is swept up and then down. 
Each recorded point is averaged over 100 ms, thus includ- 
ing 1000 excitation cycles. Such alternating excitation 
should, in principle, eliminate hysteresis in the response 
of the DC-SQUID, provided that all characteristic time 
scales of the device are much shorter than the excita- 
tion period. Note that, the thermal relaxation rate of 
a local hot-spot, created in a nano-bridge having good 
thermal coupling to the substrate, is typically on the or- 
der of gigahertz In contrary to our expectation, 
the measurement of hysteretic current-voltage character- 
istic behavior under alternating current excitation indi- 
cates that, there is a much longer characteristic thermal 
time-scale, on the order of millisecond, which affects the 
DC-SQUID dynamics. Such a long time-scale may ex- 
ist in our device because the SiaN4 membrane, to which 



the SQUID is coupled, has relatively large heat capacity, 
and also because of the poor coupling between this mem- 
brane and its surrounding, which results in a relatively 
long thermal relaxation time. 

Figure [2jb) summarizes several measurements like the 
one shown in Fig. [UJa). It plots the increasing and de- 
creasing threshold currents (i.e. the values of applied 
currents corresponding to voltage jumps in the increasing 
and decreasing current sweeps respectively) as a function 
of the lock-in amplifier carrier frequency. These thresh- 
old currents are approximately equal one another only up 
to frequencies of about 500 Hz. The hysteretic nature of 
the DC-SQUID emerges when the excitation frequency 
further increases. The decreasing threshold current falls 
to lower values while the increasing one only slightly de- 
graded. This trend continues up to frequencies of about 
40 kHz, where both threshold currents cease to depend 
on the excitation frequency. 

In order to account for the above results, we extend the 
model in Ref. Q of hysteretic, metastable DC-SQUID 
to include slow variations of the SQUID temperature. 
According to the model, the equations of motion DC- 
SQUID are the following two equations for the NBJJs 
phases 71 and 72 

7i + Pd7i + (1 + "0) V (9) sin 71 

+ J-(7i - 72 + 2tt$ x /$ ) =I x /I c0 + g n , (1) 
Plo 



72 + Pe>72 + (1 - a o) V (©) sin 72 

-J- (7i - 72 + 2tt$ x /$ ) = / x //c0 + 9n, (2) 
Plo 

where the overdot denotes a derivative with respect to 
a normalized time parameter r = uj p \t, where w p i is the 
DC-SQUID plasma frequency. In addition, I x is the bias 
current, <!> x is the external magnetic flux, $0 is flux quan- 
tum, /3d is the dimcnsionlcss damping coefficient, and Eq 
is the Josephson energy. The critical current of the DC- 
SQUID is temperature dependant, thus we employ the 
notation I c q for the SQUID critical current at the base 
temperature Tq, defined by the coolant. The dimension- 
less parameters /3lo and ao characterize the DC-SQUID 
hysteresis and asymmetry at the base temperature To, 
respectively. The dimensionless factor g n is a noise term, 
which is neglected in our numerical calculations. 

Assuming that, the temperature of the SQUID, T, 
has no spatial dependence. The term y (9) expresses 
the dependence of the NBJJs critical currents on the 
temperature. It is given by y (9) = y (9) jy (9 ) [HI, 
where 9 = T/T c , 9 = T Q /T C , and T c is the critical 
temperature of the DC-SQUID. The function y is given 

by V (9) = (1 - 9 2 ) 3/2 (1 + 9 2 ) 1/2 . Using the notation 
pc = 2itCtT c /&oI c o and fin — Ht/Ct^ p i, where Ct is 
thermal heat capacity and is heat transfer rate, the 
SQUID heat balance equation reads 
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FIG. 3: Intermittent behavior, (a) Measured average squared 
voltage of the DC-SQUID, (hi) Time traces of the DC-SQUID 
voltage taken at the working points marked by the corre- 
sponding number i in panel (a). 
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Simulation results showing the dependence of the hys- 
teretic behavior of the DC-SQUID on the excitation fre- 
quency are shown in Fig. [He). This panel shows the in- 
creasing and decreasing threshold currents as a function 
of the excitation frequency, normalized by the thermal 
relaxation time. These results show qualitative agree- 
ment with the corresponding experimental results shown 
in Fig. Hb). 

Figure [3] shows intermittent behavior of the SQUID. 
In this measurement the DC-SQUID is biased near the 
threshold of spontaneous oscillations using alternating 
current having frequency of u x /2tt = 70 Hz. The volt- 



age across the SQUID is measured as a function of time 
for various values of magnetic flux applies through the 
SQUID. Panel (a) shows the average value of the squared 
voltage versus the magnetic flux. Panels (hi) show volt- 
age time traces measured for the corresponding marked 
points in panel (a). The low amplitude voltage oscilla- 
tions are measured due to the existence of a parasitic 
serial resistance in the measurement wiring. Each time 
the alternating current drives the SQUID into the oscil- 
latory zone it responses with a voltage spike measured 
on top of the low amplitude parasitic voltage. Spikes can 
be positive or negative depending on the polarity of the 
driving current. The measured traces show intermittent 
behavior, in which the spikes are generated in bunches. 
This behavior suggests that in this region the system be- 
comes thermally bistable. The bistability is realized by 
switching between the cold locally stable state, where the 
DC-SQUID is in the superconducting state, and the hot 
locally stable state, where the SQUID becomes normal. 
Switching between these states is randomly triggered by 
external noise. 

In conclusion, we have studied nano-bridge based DC- 
SQUIDs fabricated on a S13 N4 membrane. We find that 
the response of the DC-SQUID to a low frequency al- 
ternating excitation current is hysteretic, and that the 
strength of the hysteresis depends on the excitation fre- 
quency. In addition, when the SQUID is biased near a 
threshold of spontaneous oscillations, the measured volt- 
age has an intermittent pattern, which depends on the 
applied magnetic flux through the SQUID. Such hys- 
teretic response degrades the performance of the sus- 
pended SQUID, and questions the effectiveness of using 
suspension as a method for noise reduction. 
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